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METHOD TO VERIFY THE EFFICIENCY OF ANTI- JAMMING FOR A 
COMMUNICATIONS SYSTEM 
BACKGROUND OF THE INVENTION 

The invention relates to a method to verify the efficiency of anti-jamming 
by antenna processing in one or more space communications on board a 
geostationary satellite, as well as its implementation from the ground. 

The invention can be applied in anti-jamming for space 
telecommunications and is a tool of assistance in optimizing the planning of 
connection bit rates in a theatre of operations depending on jamming conditions. 

At present, anti-jamming by antenna processing is the most efficient way 
to protect one or more space or radio communications links against hostile jamming 
units. Anti-jamming by antenna processing consists in implementing what is called 
an adaptive antenna at reception. The chief property of this adaptive antenna is that 
it matches its radiation pattern in real time to the received signal, setting pattern holes 
toward the jamming units while at the same time maintaining sufficient gain in the 
direction of the link or links to be protected as can be seen in figure 1. This result 
can be obtained from a minimum amount of information on the links to be protected 
such as knowledge of the position of the transmitters, the theatre of operations or the 
learning sequences conveyed by the transmitters without a priori knowledge of the 
jamming units present. However, in certain cases, the a priori estimation of the 
positions of the jamming units may be advantageously used by the adaptive antenna 
so as to simplify the processing operations. 

Figure 2 shows an adaptive antenna with a purely spatial structure. It is 
formed by an network of sensors Ci or radiating elements, a set of digital or analog 
reception chains CRi, downstream from the sensors, a set of adaptive filters Fi with 
one complex coefficient per filter whose role is to carry out the phase and amplitude 
weighting of the signals coming from the difference sensors before summation, and 
an adaptive algorithm A whose role is to carry out the real-time matching of the 
coefficients of the adaptive filters so as to optimize a criterion as a function of the 
information available a priori on the signals of interest and therefore the application. 
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The adaptive antenna can be implemented in an analog, digital or hybrid 
way. In the first case, the weightings are computed and applied analogically while, 
in the second case, they are computed and applied digitally. In the third case, the set 
of complex weightings is computed digitally and copied analogically before 
summation. 

For a digital implantation, the adaptive filters are formed by complex 
weighting operations whereas, for an analog implantation, these filters are formed by 
the cascade connection of a phase-shifter and a variable attenuator or a hybrid 
quadrature as can be shown in figure 3. In the context of space communications, 
when there are no jamming units, the set of weightings synthesizes a coverage (or a 
spot) on the earth, centered on a given point and having a certain surface area as can 
be seen in figure 4. In general, the coverage is characterized especially by the 3 dB 
width of the beam formed by the set of weightings. According to the size of this 3 dB 
lobe width or width of the illuminated surface area of the earth, we may speak of 
theatre, regional or global coverage, the latter corresponding to coverage of the entire 
earth. The working stations are deployed inside a coverage considered for a given 
mission and communicate together and/or with the mainland by satellite. 

A jamming operation from one or more terrestrial regions jams the useful 
uplinks (from the earth to a satellite) and it is the role of the adaptive antenna 
precisely to carry out anti-jamming on the links by creating antenna pattern holes 
toward the jamming units, located outside or within the coverage and picked up 
respectively by the minor or major lobes of the antenna. 
SUMMARY OF THE INVENTION 

The invention relates to a method for the verification of the efficiency of 
the anti-jamming, by adaptive antenna, of the uplink of one or more space 
communications links as well as its implementation from the ground. 

The invention relates to a method of anti-jamming in a communications 
system comprising several sensors or adaptive antennas. It is characterized by the 
fact that it comprises at least the following steps: 

• estimating the mean power of the output of the communications system, 
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• estimating the respective power values Pu or P'u, of a station u, the antenna noise 
Pa or P'a, the thermal noise PT, or P'T, 

• estimating at least one of the following ratios: 

P v 
J«JS<e^= (g; ; P p ) / (g; ; 

— m 



Jtot/Stot = ( g W <Z P " > (22) 

p=l «=1 



with p = the jamming unit 

= sum of the power values of the residual sum/sum of the power values of the 
stations on the reception band B. 

JterlS ^- (g; ; (23) 

p 

JtotlSu = ( g P p)I P« (23) 
p=i 



= sum of the power values of the residual jamming units/power of the station u 
in the reception band B. 



p 




p 

J U IS U = ( g Ppu)/P u (24) 



With Ppu = power of the jamming unit p in the reception band Bu. 
• comparing at least one of the three ratios with a threshold value. 
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The invention also relates to a system for the verification of anti-jamming 
in a communications system comprising several sensors or adaptive antennas, and a 
piloting device on the ground. It is characterized by the fact that it comprises at least 
the following elements: for a verification by channel, from the ground and for a 
5 reception band B, a computer integrated into the piloting device and an onboard 
computer, the two computers being programmed to execute the following steps: 
Communications Channel Power Measurement: Onboard function parametrized 
from the ground by the Onboard Param VAA function, 
VAA GAIN: Ground function Sol, 
10 Communications channel power measurement : onboard function, 
VAA Processing : Ground function. 

According to another alternative embodiment, the invention also relates 
to a system for the verification of anti-jamming in a communications system 
comprising several sensors or adaptive antennas, a piloting device on the ground 
1 5 comprising at least the following elements : 

For a verification by stations, an onboard computer and a ground 
computer, the computers being programmed to execute the following functions : 
Communications Channel Power Measurement : onboard function parametrized 
from the ground by the Onboard Param VAA function, 
20 VAA Gain : ground function, 

Acquisition of Communications Channel : onboard function parametrized from the 
ground by the function Onboard Param VAA, 
VAA Processing : ground function. 

The invention can be applied for example in space communications 

25 systems. 

With the proposed method, it is possible at all times to know whether or 
not the anti-jamming applied is effective. If it is not effective, the information 
coming from the method modifies the anti-jamming characteristics (choice of the 
number of type of auxiliary channels in the case of an OLS (Minor Lobe Opposition) 
30 type of processing, alternative parametrization of a pre-synthesis of zeros (PRS) 
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when a piece of a priori information on the position of the jamming units is available 
etc) to increase its efficiency. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the invention shall appear more clearly 
from the following description given by way of an illustration that in no way restricts 
the scope of the invention, along with the appended figures of which : 
Figure 1 shows a radiation pattern of the antenna after anti-jamming, 
Figure 2 is a functional diagram of a spatial-structure adaptive antenna, 
Figure 3 is a purely spatial adaptive filter for an analog implementation, 
Figure 4 shows the coverage demarcated by the beam associated with the set of 
weightings when there are no jamming units, 

Figure 5 shows a structure of the adaptive antenna for a digital implementation of the 
filters, 

Figure 6 shows a structure of the adaptive antenna for an analog implementation of 
the filters, 

Figure 7 is a functional diagram of the sequencing of the operations for the 
implementation of the system for verifying the efficiency of the anti-jamming. 
MORE DETAILED DESCRIPTION 

The method according to the invention uses especially the information, 
assumed to be available a priori, on the position and the Equivalent Radiated 
Isotropic Power (ERIP) sent out by the working stations working within the 
coverage, also called theatre information. It furthermore makes use of the 
characteristics of the active antenna used onboard the satellite and especially 
knowledge of the positions and responses of the RE (Radiating Elements) for each 
direction of space and each polarization of the incident field, the set of weightings 
used for the anti-jamming, the gain and the equivalent noise temperatures of the 
analog or a digital reception chains downstream from the sensors and, for an analog 
or hybrid layout of the set of weightings, that of the digitization chain if any at output 
of the antenna. 

Before explaining the method according to the invention, a few 
reminders are given on the signals in a communications system with anti-jamming. 
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A. Signals at output of a communications BFN (Beamforming Network) 

It is assumed that each of the TV sensors Ci of the network or array of the 
figure 2 receives the contribution from U useful sources, coming from the theatre of 
operations, from P jamming units disturbing the communications and from a 
background noise. It is assumed that all the signals are narrowband signals for the 
network or array of sensors. 

Al. Expression for a digital implementation of the adaptive antenna 

Figure 5 shows the structure of the adaptive antenna in the case of a 
digital implementation of the adaptive filters. 

The N sensors of the network correspond either to REs, or to sub- 
networks or sub-arrays preformed in analog mode. In the context of a digital 
implementation, the vector, x(t), of the envelopes of the signals brought to the point 
PI of the figure 5 is written as follows at the instant t 

v P 




u P 



where b a (t) is the noise vector at the point PI coming from the network or array of 
sensors or antennas (external noise + thermal noise of the RF reception chains), bj(t) 
is the thermal noise vector of the digital chains brought to ?l,j p (f) and J p correspond 
respectively to the complex envelope and to the direction vector of the jamming unit 
p, sjf) and S u respectively correspond to the complex envelope and to the direction 
vector of the station u. 

In the general case pertaining to any unspecified sensors, the component 
n of the direction vector S u is given by 

Sun = m,tl«)exp(-j *„#•„) (2) 
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where k u and r|„ are respectively the wave vector and the polarization parameters of 
the station u, r n is the position vector of the sensor n anAf n (k u , n„) is the complex 
response of the sensor n in the direction k u for the polarization r\ u . 

On the above assumptions, the complex envelope at the instant nT e , y(n), 
of the sampled output of the jamming-protected communications BFN associated 
with the set of weighting operations w, is written as follows : 

><n) ~; A w ^-G^-m-^ (3) 

v p 

u~ 1 p- 1 

y(n) & yJ G num x(n) = (3) 
U P 

s u (n) w< G mm S u +Yj j p (n) w t G num J p + w^G mm b a (n) + w*G num b^n) 



where T e is the sampling period, G mm is the diagonal matrix (NxN) whose diagonal 
elements are the gains of the digitization chains. 
20 A2. Expression for an analog or hybrid implementation of the adaptive antenna 

Figure 6 shows the structure of the adaptive antenna in the case of an 
analog or hybrid implementation of the adaptive antenna, namely for an analog 
application of the adapter filters. 

The N sensors Ci of the corresponding networks are either REs or sub- 
25 networks preformed in analog mode. In the context of an analog implementation, the 
vector, x(t), of the envelopes of the signals brought to the point PI of the figure 6 is 
written as follows at the instant t 
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*(0 = 2] Su(t)S u + £ Jp{t)J p + 



(4) 



where b a (t) is the noise vector at the point PI coming from the network of active 
sensors (external noise + thermal noise of the RF reception chains) and/or the other 
parameters defined in the above paragraph. 

On the above assumptions, the complex envelope, y(n), of the sampled 
output of the jamming-protected communications BFN associated with the set of 
weightings w, is written as follows : 



15 bmi 



y(n) - a {w* G x{n) + bj{n) } 

U P 

20 = a { 2] s u (n) ^G^+J] J p (n) w^GJ p + w^G b a (n) + bj(n) } (5) 

u=l P =i 



where G is the diagonal matrix (NxN) whose diagonal elements are the gains of the 
25 RV chains, w is the vector of the analog weightings, a the gain of the digitization 

chain of the output of the BFN and bj(n) is the thermal noise of the digitization chain 

of the output brought to the point P3. 

In practice, the matrix G is generally known for a reference temperature 

T 0 and is denoted G 0 . For an antenna temperature, T Ant , the matrix G is no longer 
30 equal to G 0 but takes the value 
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G = [G Q Z + (Tji nt — T 0 ) 8G 2 1] 



(6) 



where 8G is a coefficient of variation of the gain in amplitude of the RF chains with 
the temperature and I is the identity matrix. 
5 B. Power of the output of a jamming-protected communications BFN 
Bl. Expression for a digital implementation of the adaptive antenna 

Assuming that all the signals are decorrelated from one another, from the 
equation (3), we deduce the power of the output of the communications BFN in the 
case of a digital application of the adapter filters, given by 



zr^ — 

*-=■ 4 



*y ~ < E[| yinf] > = w t G num R x G n J w (7) 
u P 

Z *u\yJ G nwn s u \ 2 +Y, M ^ G num J p | 2 + (n a + nx) ^G num G n J w 



where O corresponds to the operation of temporal averaging on an infinite horizon 
of observation, ^ ~ A <£[*(/*) x(nf ]>- ^-<EW^) f ]> is the averaged matrix 
of correlation of jc(»),^ w -; A <E[|j w fo)f }>- - <E[|^(«)| 2 ]> j s me mean p0W er 
of the station u picked up by an omnidirectional KE,-n p ~; A <E[|fc fo)f ]> 

- <EH/>Cw)| .]> i s the mean power of the jamming unit p picked up by an 
omnidirectional RE, n a and n T , such that <E[b a (n) b a (rrf]> = rj a I and <E[bj(n) 
- tit I, are the equivalent mean power values per sensor brought to the point 
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PI in terms of antenna noise and thermal noise respectively, assumed to be spatially 
white. 

In introducing the power values P u , P p , P a , P T , respectively of the station 
u, the jamming unit p, the noise of the antenna a and the thermal noise of the 
5 digitization chains at output of the communications BFN, respectively defined by : 



P u = n u \vJ G„ um S u \ 2 
P P = K P \>J G num j p ? 

Pa = Tla W^G num G n J W 
P T = tit W*G num G nu J W 

the expression takes the following form : 



(8) 
(9) 
(10) 
(11) 



— m 

15 

U P 
•*9 = 2 P u + J] Pp + P* + P T 



20 The power values r| a and r| T are given by 

r\ a =kT a B (13) 
r] T =kT T B (14) 
where k is the Boltzman's constant, B is the reception band and T a and T T are the 
temperatures of the equivalent antenna noise and thermal noise per sensor at PI. The 
25 equivalent thermal noise temperature at PI per sensor, T T , is computed from the 
ambient temperature, T amb , and from the noise factors of the elements of the 
digitization chain for the sensor considered. In practice, the equivalent antenna noise 
temperature at PI is generally known for a reference temperature T 0 and is denoted 
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T a0 . For an antenna temperature, T Ant , the noise temperature T a is no longer equal to 
T a o but takes the value 



T a = T a0 + (T Ant -T Q )8T (15) 
where 8ris a noise temperature gradient relative to the temperature of the antenna, 
known a priori. 

Furthermore, the power n u of the station u is linked to its ERIP, ERIP(w), 
by the following expression 

7t M = ERIP( M ) ( X I4ur u f (16) 
where X is the wavelength of the carrier wave, and r u is the distance between the 
station u and the satellite. A similar relation links the power n p of the jamming unit p 
and its ERIP, ERIP(p). 

B2. Expression for an analog implementation of the adaptive filters 

Again assuming signals that are decorrelated from each other, from the 
expression (5), we deduce the power of the output of the commumcations BFN in the 
case of an analog application of the adaptive filters expressing 

sr — ^ < PH v ( Vl i a i > - I rY l^ f w*C li C^ w 4- n -r i 

y~ j ^ [Da ) I J ' \ \ \ vj ~ — =?= — rj j -j- 

U P 



n y & <E[|j(")| 2 ]> = N 2 {w f GR x G*w + r\ T ) 



= N 2 {X ^l^^^| 2 + Yi "p\^ gj p\ 2+ i\*JGd*w + t!t} (17) 



where n a , such that <E[A a (/i) 6 a (n) t ]> = r| a I, is the mean power, at the point PI, of 
noise per sensor coming from the active network (external noise + thermal noise of 
the reception chains), ^ ~; A <E[\b^ n)f^>- *it - <E[\bj(n)\ 2 ]> j s t h e mea n power 
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of thermal noise coming from the digitization chain brought to P3. The quantities r\ a 
and tit are defined respectively by (13) and (14) where T a is the equivalent thermal 
noise temperature per sensor of the active antenna at PI and where T T is the 
equivalent thermal noise temperature coming from the digitization chain and brought 
to P3. Similarly, the power values % u and % p are related to the ERIP by the expression 
(16). 

In introducing the power values P' u , P' p , P' a , P' T , respectively of the 
station u, the jamming unit p, the noise of the antenna and the thermal noise of the 
digitization chain at output of the communications BFN defined respectively by : 

P' u = \a\ 2 n u \w f GS U \ 2 (18) 
P' p = \a\\ p \w^GJ p \ 2 (19) 
P' a = |a| 2 n a w*G G f w (20) 
P' T = |a| 2 r| T (21) 
the expression (17) takes the form (12). 
Principle of the invention 

The steps of the method according to the invention rely especially on the 
following idea : from an estimation of the mean power, % y , the output of the 
communications BFN and the estimates of the quantities P u , P a et P T , P' u , P' a and 
P't , the method makes it possible to estimate the efficiency of the anti-jamming. 
This is done especially by estimating different residual jamming unit/stations ratios, 
in the reception band, at output of the communications BFN. 

For example, according to an exemplary implementation, the method 
uses three residual jamming unit/station ratios whose values make it possible to 
evaluate the efficiency of the anti-jamming or of the set of weightings w considered 
at output of the communications BFN. The three ratios considered here below in the 
document correspond to : 

- the ratio of the power values respectively of the sum of the residual 
jamming units to the sum of the stations in the reception band B, 
hereinafter called (J/S per channel) and referenced J t0 t IS tot 
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- the ratio of the power values respectively of the sum of the residual 
jamming units to the power of the station u in the reception band B, 
hereinafter called (J/S u per channel) and referenced J tot IS u 

- the ratio of the power values respectively of the sum of the residual 
jamming units to the power of a station, in the band B u of the station, 
hereinafter called (J/S„ per station or per link) and referenced, for the 
station u, J/S u . 

These quantities are defined respectively by : 



(g; ; P p )l (g; ; P u -+ 




-(33) 



p u 



(22) 



P-l «=1 



P 



( gj 



p 



JtotISu = ( ^ Pp)/ P„ 



(23) 



-(24) 



P^i 



p 



JulS u = (^ P PU )/P U 



(24) 



where P p „ is the power of the jamming unit p in the band B u . 

The quantity P u , P' u respectively defined by (8) or (18), is estimated 
from a priori knowledge of the theater of operations (ERIP and position of the 
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working stations), the center frequency of the reception band, the responses of the 
sensors of the network, the set of weightings w as well as the gains of the reception 
and digitization chains, G num , G, a, known on an a priori basis or computed by (6) 
from the temperature of the antenna and of the parameter 8G. 

The quantity P a , P* a respectively defined by (10) or (20) is estimated 
from a priori knowledge of the set of weightings w, the gains of the reception and/or 
digitization chains, G num , G, a, known on an a priori basis or computed by (6) from 
the temperature of the antenna and the parameter 8G (&G is a coefficient of variation 
of the gain in amplitude of the RF chains with the temperature) as well as from the 
equivalent noise temperature of the antenna, T a , in PI (itself a function of the 
temperature of the antenna, T ant ), the reference temperature T 0 , the temperature of the 
antenna noise, T aQ at PI at the temperature T 0 and the variation in noise temperature, 
8r, with the temperature. 

Finally, the quantity P T , P' T defined by (11) or (21), is estimated, for a 
digital implantation, from a priori knowledge of the set of weightings w, gains of the 
digitization chains, G num , as well as the temperature T r of the equivalent thermal 
noise per sensor in PI. For an analog implantation of the filters, the quantity P T is 
estimated from the knowledge of the gain, a, the digitization chain at output of the 
BFN and the thermal noise temperature of this chain brought to P3, T r . In both 
cases, the quantity T r is estimated from the ambient temperature T amb and from the 
noise factors of the elements constituting the digitization chain or chains. 
C. Estimation of the J/S at output of the communications BFN 

For a digital implantation of the filters, the estimation of the ratios 
defined by the expressions (22) to (24) necessitates the estimation of the quantities 
*y, p u, P a and P T defined respectively by (7), (8), (10) and (11) and, for an analog 
implantation of the filters, it necessitates the estimation of the quantities P' u , P' a and 
P\ defined respectively by (17), (18), (20) and (21). 
CI. Estimation of % y 

The method estimates the mean power n y of the output of the 
communications BFN from a number K of samples, y(k), 1 < k < K, of this output. 
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For a sufficient oversampling factor, an asymptotically unbiased estimator of this 
mean power is given by : 



, — 2j ; f 



(25) 



10 This estimator becomes consistent for stationary and ergodic outputs and 

cyclostationary and cycloergodic outputs. 
C2. Estimation of P u 

The method estimates the values p£ H , P;~' w - ■ P u > P\ of the 

power P u , P' u defined by (8) or (18) in using firstly the a priori knowledge of the 
1 5 parameters w and G num for a digital application of the adaptive filters and |a| 2 , w and 
G for an analog application of these filters and, secondly, the estimation of the 
parameters % u and S u . 

The applied set of weightings w is permanently known while the matrix 
gain G num and scalar gain |oc| 2 of the digitization chains are parameters adjustable 
20 from the ground by the operator so as to optimize the use of the dynamic range of the 
ADC or ADCs as a function of the jamming environment. The matrix G of the gains 
in amplitude of the analog reception chains is mastered, through the expression (6), 
from the knowledge of the matrix G 0 of the gains for the reference temperature T 0 , 
the parameter 8G and the permanent control of the temperature of the antenna T Ant . 
25 The mean power, n u , of the station u, received by an omnidirectional 

sensor is estimated by the expression (16) where the ERIP, ERIP(w), of the station u 
is known a priori and listed in a mission plan, where X is deduced from the 
frequency channel and where r u is deduced, for a geostationary satellite, from the 
position of the station u on the earth. 
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Finally, the direction vector S u , whose component n is given by (2), can 
be deduced from a priori knowledge of the positions, r„ of the sensors of the 
network, of the wave vector k u through the position of the station u, the polarization, 
r|„ of the station u and the complex responses f n (k u , r\ u ) of the sensors for the wave 
vector k u and the polarization r|„. 
C3. Estimation of P a 

The method estimates the values P^ aT Pp- s - p * p, a of the 

power values P a , P' a defined by (10) or (20), in using firstly the a priori knowledge 
of the parameters w and G num for a digital application of the adaptive filters and |a| 2 , 
w and G for an analog application of these filters and, secondly, the estimation of the 
parameter r| a . 

The mastery of the parameters w, G num , G and |a| 2 is discussed in the 
previous paragraph. The estimation of the power, n a , of the noise of the antenna per 
sensor at the point PI is computed by the expression (13) where the equivalent noise 
temperature of the antenna, T a , at PI is obtained by the expression (15) from the a 
priori knowledge of the reference temperature T 0 , the antenna noise temperature, T a0 , 
at PI at the temperature T 0 , the variation of the noise temperature, ST with the 
temperature and the permanent measurement of the temperature of the antenna T Ant . 
C4. Estimation of P T 

The method estimates the values P^^-PiP-^- Pt,P't' of the 

power P T , P' T defined by (1 1) or (21), and requires, firstly, a priori knowledge of the 
parameters w and G num for a digital application of the adaptive filters and |cc| 2 for an 
analog application of these filters and, secondly, the estimation of the parameter n T . 

The control of the parameters w, G mm and |oc| 2 is given in the paragraph 
B. The power, n T , is estimated from the expression (14) where T T is the temperature 
of the equivalent thermal noise of a sensor digitization chain brought to PI, for an 
application of the adaptive filters in digital mode, and of the digitization chain of the 
output of the BFN brought to P3, for an application of the filters in analog mode. In 
both cases, the quantity T r is estimated from the ambient temperature T amb and from 
the noise factors of the elements constituting the digitization chain or chains. 
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With the different values n y , P u , P a and P T having been estimated, 
according to the method, at least one of the three vatios^te^^te^^t^-^ur^F 
J tot I Stot , J m I Su, JISu is estimated. The expressions have been 

given for the case of a digital application of the adaptive filters and remain valid in 
exchanging the letters, Pu, P a and P T , by the letters P «, P' a and P' T } for an analog 
application of the adaptive filters. 
C5. Estimation of J tot /S tot 

From the above estimations, we deduce an estimation,-^ ^ / 
I™ 1 §ob of the ratio J tot IS tot defined by (22), given by 



Pr-»-J (26) 

u u 
Jtot/Stot = (% - 2] h - Pa ~ Pt )/ (^ ) (26) 



C6. Estimation of J tot IS U 

From the above estimations, the method deduces an estimation,-^,^ 
§Jz hT . J tot I S u , of the ratio j t j Su defined by (23), given by 

u 

£; ; — p£»-^*—p£rV-p£u 



c/ 

= (^ - J] Pu - Pa~ Pt)/ h (27) 



-(37) 
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C7. Estimation of JIS U 

The estimation, jf/S?^ J/s «' of the ratio J /S u defined by the 
expression (24), necessitates the estimation of the total power of residual jamming 
units in the band B u of the working station u. This estimation necessitates the 
following operations : 

- reception of the samples, y(k), of the output y(t) of the communications 
BFN, 

- bandpass filtering of the samples around the band B u . The samples y u (k) are 
obtained, 

- estimation of the power of the output filtered by (25) where y u (k) replaces 
y{k). We obtain*^ 

estimation of the power values of antenna noise and thermal noise 
respectively at output of the BFN in the band B u . These quantities are 
computed, from the equivalent noise temperatures computed here above, by 
the expressions (13) and (14) respectively, where B is replaced by B u . We 
thus obtain Pf 1 ^ and P;~q ^ PgumdPn,, 

- Computation of the power of the stations v other than the station u in the 
band B u at output of the communications BFN. The approach is that of the 
step B but one in which, for each station v different from u, the ERIP used 
in the computation of-*Ar- 2*- is that of the station v in the band B u . 
Thus the quantities PP m — P™ are obtained, 

- computation of the ratio J-~/S;~ u — JIS U , ' by the expression : 



Ji/sr t 



-£tr 



v*u 
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With the estimate of at least one of the three ratios being known, the 
method compares the estimated value or values with a threshold value Vs. 

If the value found is above this threshold value then, according to the 
invention, a message is sent on the inefficiency of the anti-jamming. If not, the 
message informs, for example, an operator that the jamming efficiency is sufficient. 

The threshold values take account firstly of the permissible jamming 
power per station or per channel to carry out the demodulation of the stations and, 
secondly, the precision of estimation of the previous ratios. The computations of 
precision made in the paragraphs Di show that, for jamming unit/station ratios at 
output greater than 0 dB, the precision of estimation of these ratios by the proposed 
method is very high whereas this precision decreases with the ratios between the 
jamming units and the output signal. In this context, it may be considered that the 
anti-jamming is not efficient if the ratios between the jamming units and the output 
station exceed 0 dB. 

Thus, the control of the precision with which the estimators given here 
above estimate the different Jamming unit/Station rations considered enables 
especially an efficient operational exploitation of these estimators. For this reason, 
the method may comprises a step for determining the precision of each of these three 
estimators. 

Dl. Precision of estimation of % y 

The estimate-^f^- Ky_ is related to % y by the following expression : 




Tty - %y(l+ Ally ) 



(29) 



where An y characterizes the error on the estimation of A^. the previous expression 
expressed in dB becomes 



foy)dB = (ttj,)dB + lOlogio(l-HA^) & (jtj,) dB + H^y)dB 



(30) 
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Assuming that the samples y(k) are independent (all the sources are 
spread out in the reception band B), stationary, Gaussian, the estimator (25) is not 
biased (E[7t;^ =^- ( E W = and has a variance : 

Vu&£&=-*?-l-K (34) 

Var[^] = 7t//Z (3!) 

i.e. a mean standard deviation of 

^bg&^+fe (32) 

gfe] = KyHK (32) 

Thus, in 99% of cases, the estimator-*^- .% is such that 

- %(1 3/ViQ < ti;~^ < ^(1 +3/V^) (33) 

7^(1 -3 < TCy < 7Ty(l-f3/Vg) (33) 

where Arc^, is a centered random variable, that is quasi-Gaussian for K as 
a great value and having a mean standard deviation 1 Nk. Thus, in 99% of the cases, 

-3/<K < An y < 3 / V/sT (34) 
101ogio(l -3 / VisQ < A(7t J ) dB < 101ogio(l + 3 / V#) (35) 

Digital application: 
For K= 1000, we obtain 

- 0.4 dB < A(% y ) dB < 0.4 dB i.e. precision of ± 0.4 dB. 



D2. Precision of estimation of P u 
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From the expressions (8), (16) and (18), we deduce the expression of the 
power P u at output of the communications BFN for an application of the digital and 
analog filters respectively, given respectively by 



P u =WP{u){X/Anr u ) 2 \^G num S u \ 2 
P'« = ERIP(w) ( X 1 4 n r u f \a\ 2 \>J GS U | 2 



(36) 
(37) 



This means that the power of the station u,-P-t- ur p » reconstructed 
from the information on the mission, can be written as follows for a digital and 
1 0 analog application, respectively, of the filters 



_ A A 

fit:. — " P. (\ + ,\P \ — P I 



P u ^ PuQ- + AP U ) Z g«(j + AERIP(»)) (1 + A| yJ G num S u \ 2 ) 



(38) 



J, w (l + AP' U ) 



K Z P '«( l + ^jj = gjfO + AERTPQ/)) (1 + A|q| 2 ) (1 + A| G S u j 2 ) (39) 



where the quantities AERIP(w), A|a| 2 A^G,^/ et A^GS*/ are the errors in the 
knowledge, respectively, of ERIP(u), |cc| 2 , | w*G num S u | 2 et | w*G S u | 2 . 

From the previous expressions, we deduce that of-p£ — p in dB, 

25 given by 

Wide— M^dB-^44k^ 

(PJdB = (P u )dB + 101ogio(l+AP„) * (PJm +A(P u ) dB (40) 



30 
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where, for an application of the filters in digital mode 
A(P„)dB = A(ERIP( W )) dB + ^G num S u f) AB 



(41) 



whereas, for an application of the filters in analog mode 
5 A(i"„) dB = A(ERIP( M )) dB + A(|a| 2 ) dB + A(| w^GS u | 2 ) dB (42) 

Dieital Application : 

For example, if it is assumed that 

- the uncertainty on the ERIP of the stations, A(ERIP(u)) dB , is ±2 dB, 

10 -the uncertainty on the gain of the digitization chain of the output of the 

BFN, A(/ap) dB , is ± 0.5 dB (because of the drifts in temperature and the effective 
application of the gain), 

- the uncertainty, A(| w*G num S u | 2 ) dB , is ±1 dB because of± 0.5 dB of uncertainty 
on the gains of the digitization chains and ± 0.5 dB of uncertainty on the 

1 5 components of the direction vector S u because of the uncertainties on the position 

of the station and on the responses of the sensors, 

- the uncertainty, A(| w^G S u | 2 ) dB , is+ldB for the same reasons as above. 

The value A(PJ dB = ± 3 dB is obtained for a digital application of the filters and 
^(PJdB = ±3.5 dBfor an analog application of the filters. 
20 D3. Precision of estimation of P a 

From the expressions (10), (13) and (20), we deduce the expression of 
the power P a at output of the communications BFN for an application of the digital 
and analog filters respectively, given respectively by: 



from the information on the reception chains, the set of complex weightings and the 
antenna noise temperature can be written as follows for a digital application and an 
analog application, respectively, of the filters : 



25 




This means that the antenna noise power,^£^ a — P a > reconstructed 
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(45) 



-(45) 
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n^ G G^ w) ) 

P'a ~ P'S + AP' a ) = P' a (l + &T a ) (1 + Alaft (1 + A(w^G G t w)) (46) 



pertaining to the knowledge, respectively, of T a , |cc| 2 , w t G num G nu J w 
and w^G G^ w. 

From the above expressions, we deduce that of-^f^— Z*. m dB» given by 

&)«m = (^)dB + 101ogio(l+AP a ) £ (P a )dB + A(^) aB (47) 

where, for an application of the filters in digital mode 

A(P a ) d B = A(T a ) dB + A(w^G num G n J w) dB (48) 

whereas for an application of the filters in analog mode 

MP'a)dB =A(r a ) dB + A^GG^jds + A(|a| 2 ) dB (49) 

Digital avvlication : 

For example, if it is assumed that : 

- the uncertainty on the antenna temperature is ±0.5 dB, 

- the uncertainty on the gain of the digitization chain of the output of the 
BFN, A(/af)4 B , is ±0.5 dB (because of the drifts in temperature and the effective 
application of the gain), 

- the uncertainties, A(w^G num G m J w) dB and A(w T G G t w) dB are ±0.5 dB because 
of± 0.5 dB of uncertainty on the gains of the RF and digitization chains. 

A(Pa)dB = ±1 dBis obtained for a digital application of the filters and A(PJd B = ± 1.5 
dB for an analog application of the filters. 
D4. Precision of estimation of P T 
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From the expressions (11), (14) and (21), we deduce the expressions of 
the power P T at output of the communications BFN for an application of the digital 
and analog filters respectively, given respectively by 

P T = kT T B w*G num G n Jw (50) 
P* T = kT T B |cc| 2 (51) 

where T T has a different sense depending on the nature of the implementation. This 
means that the antenna noise power,-^^ reconstructed from the information 

on the reception chains, the set of complex weightings and the thermal noise 
temperature at PI for a digital implantation and at P3 for an analog implantation is 
written as follows for a digital application and an analog application, respectively, of 
the filters : 

P T & Pj(l+ AP T ) = Pjjl + A7» (1 + A(w^G num Gnum wj) (52) 

— m 

P' r ^ P'jjl +AP' T ) = P'i(l + A7»(1 + Ajaj 2 ) (53) 



where the quantities AT T , A|a| 2 and A(w t G WKW G nu J w) are the errors relating to the 

knowledge of T T , \a\ 2 and w^G num G nu J w respectively. 

From the above expressions, we deduce that of p£$- Pj_ in dB, given by 
(^^m-=-m m + lO\og^(^^^-(P^ m -^^^ m (§4) 

A)dB = (Pr)dB + 101ogio(l+AiV) =, A (Pr)dB +A(Pr)dB (54) 

where, for an application of the filters in digital mode 

A(P r ) dB = A(r r ) dB + A(w^G num G n J w) dB (55) 
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whereas, for an application of the filters in analog mode 

A(PV)dB =A(7» dB + A(|a| 2 )dB (56) 
Digital application : 
For example, if it is assumed that : 

- the uncertainty on the antenna temperature is ±0.5 dB, 

- the uncertainty on the gain of the digitization chain of the output of the 
BFN, Af/afjds , is ± 0.5 dB (because of the drifts in temperature and the effective 
application of the gain), 

- the uncertainties, A^G^ G nu J w) dB are + 0.5 dB because of ± 0.5 dB of 
uncertainty on the gains of the RF and digitization chains. 

A(Pa)dB = ±1 dB is obtained for a digital as well as an analog application of the filters. 
E. Precision of estimation of S tot 

The estimation^f^de-^ -SW.de 5^ j s wr itten as follows : 

U ; u 
~; A SO^^Sm) ~ ; ~ ; P U &±^P„) 

(§7) 

u u 



giving in dB 

(StoddB = (S tot ) dB + lOlogioCl+A^) = A (S tot ) d B + HS,oddB (58) 

where 

v 
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u 

AOWdB = 101og 10 (l+ - ) (59) 

Z p» 

5 Digital application : 

For example, if it is assumed that the precision on the power of the stations is identical 
for all the stations, MStodtiB ~ A(P„) dB » ±3dB or ± 3.5 dB is obtained according to the 
nature of the implantation.. 

1 0 F. Precision of estimation of J tot 

The estimations^- J joh_ 0 fJ tot is written as follows 



, (1 + Atc^ + AP U ) + AP a ) J»j <1 + AP^ 



— m 

that is, in dB 

— m 

= Meg_ 
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u 

%ot - Jtoti} +A/ tof ) = % y - J] h-Pa-Pr (60) 

(7 

= ny (1 + Arc, ) - ^ A(l + AP„) - i^l + AP fl ) - Pj(l + AP T ) (61) 

u=\ 

that is, in dB 

OtodtB = (4tf)dB .+ 101ogio(l +AJ to/ ) ^ (J to ,) dB + A(J to ,) dB (62) 
where 

1/ 

A(J tot ) dB = 101og 10 (l+ > (63) 



From this result, it is deduced that the precision of estimation of J tot 
depends on the relative signal and jamming unit contributions at sampled output of 
the communications BFN. 

More specifically, for jamming residues that are very high before the 
stations (either because of an absence of anti-jamming or because of low- 
performance anti-jamming, when there is high-level jamming at input) it is deduced 
from (63) that the precision on J tot is close to the precision on % y . 

Digital application : 

In these conditions, A(J tot ) dB « A(n y ) dB « + 0.4 dB. 

By contrast, for jamming residues that are very low before the stations 
(either because of an absence of jamming or because of high-performance anti- 
jamming) the total power is close to that of the working stations and the error may 
become very great. 



G . Precision of estimation of J m I S u 

The estimation, Est[J to , / S u ], of J tot I S u is written as follows 



MARKED-UP SPECIFICATION 

28 



wh e nc e (65) 

& + 1 SuMV (4 * A P u ) JO + 

that is, in dB 

^^tert-Sibfjg— — — 

W^g^ (B + A P u )] fl * 

ErtMpr/*] ^ (^/^d+A^/^)) = Xt/S u (64a) 

(1 + Any) - ]T p ui i + _ Pa(l + _ + ^ 

(64b) 

^d+AP a ) 

whence (65) 

V+WVHV (l + AP„) ]0 + "^-^--^-^ 

^ - Z ^« - ^ a - p T 

that is, in dB 

WtotlS u ). R = 
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From this result, it is deduced that the precision of estimation of J tot I S u 
depends on the relative signal and jamming unit contributions at sampled output of 
the communications BFN. 

More specifically, for jamming residues that are very high before the 
5 stations (either because of an absence of anti-jamming or because of low- 
performance anti-jamming, when there is high-level jamming at input) it is deduced 
from (66) that the precision on J tot l S u is given by : 

A(J tot /S u ) dB = I0log 10 ((l+A % )/ (1+APJ) = A(^) dB - A(PJ dB (67) 

10 

Digital application : 

In these conditions, A(J tot I S^)^ « ± 3.4 dB. 

By contrast, for jamming residues that are very low before the stations 
(either because of an absence of jamming or because of high-performance anti- 
15 jamming) the total power is close to that of the working stations and the error may 
become very great. 

H . Precision of estimation of J tot l S tot 

The estimation Est[/ to , / S tot ], de J tot I S tot can be written 

y 

Est[J {ef -/S {0 (\ ==7^ (J t0 f^S t0t ) (1 + Atftef/Sfe^) — — Jj-fefJ—y^^, ', Pj- u 

4- 

20 (6?a) 



S (67b) 

ice 

(1 + A(J t 0r /S t0 f)) ~ [g; ; P u I gj ; Pj ^APJ^- 

( 1 

-H68) 
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that is, in dB 

Wleg^ ( ( 1 + 



Est[J tot /S tot ] ^ (Jtot/Stodil+AVtot/Stot)) = Xt'Yi ^« (67a) 
Ttj, (1 + ATtj,) - ^ ^ u (l + AP„) - P a {\ + APJ - Pj(l + AP T ) 



(67b) 



whence 

(1 + A(7, 0/ / 5^)) = [ £ P u / Z /•„(!+ AP U ) ] x 



t^Ati^- J] P u AP u - P a AP a - P r AP T 
(1+ )( 



u u 

WtotlSto,)^ = I0log 10 ([£ i>„/£ P„(1+APJ]) 
U 

n y An y - Pu/^u ~ P<AP a ~ Pt&Pt 
+ I0log 10 ( ( 1 + ) ) (69) 

*y~ Z Pu-Pa-Pr 



From this result, it is deduced that the precision of estimation of J tot I S tot 
depends on the relative signal and jamming unit contributions at sampled output of 
the communications BFN. 
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More specifically, for jamming residues that are very high before the 
stations (either because of an absence of anti-jamming or because of low- 
performance anti-jamming, when there is high-level jamming at input) it is deduced 
from (69) that the precision on J tot I S tot is given by 

V V 

m 

u u 

MJtot I SW) dB - lOlog 10 ( ( 1 + A*,) [ X P u I £ P„(l + A/>J ] ) = A(A7ij,) dB + 
u u 

i° lo 8 10 ([S Pu'T. ^(1 + AP U )]) (70) 



which gives (67) if all the stations have the same precision. 
Digital application : 

In these conditions, from the above examples, A(J tot I S to t) dB « ± 3.4 dB. 

By contrast, for jamming residues that are very low before the stations 
(either because of an absence of jamming or because of high-performance anti- 
jamming) the total power is close to that of the working stations and the error may 
become very great. 

Exemplary implementation of the method in a communications system 

The method whose steps have been explained here above is, for example, 
used in the system comprising a base located on the ground and including a computer 
program to implement the functions described in detail here below, the base being 
linked by means known to those skilled in the art with one or more satellites 
equipped with chains such as those described in figures 5 and 6. 

Figure 7 is a block diagram of an exemplary sequencing of operations. Two 
cases of operation sequencing are possible depending on whether it is sought to 
estimate the quantities J tot I S tot and J tot I S u relative to the reception band B or, on 
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the contrary, the quantity J u I S u relative to the band of the station u. In the former 
case, the term "verification by channel" will be used and in the latter case the term 
used will be "verification by station ". 

For a verification by channel, the method carries out the operations depicted by 
a solid line in figure 7. Starting from the ground and for a reception band B, it 
executes the following functions : 

- Communications Channel Power Measurement : whose aim is to estimate 
the total power available at output of the digitization chain of the 
communications BFN for the applied set of weighting operations and send 
it to the ground. This function is an onboard function parametrized from 
the ground by the function Onboard Param VAA implemented for example 
in a computer, 

- VAA GAIN: whose aim is the optimizing, from the results of the 
Communications Channel Power Measurement function, of the gain of the 
digitization chain of the output of the communications BFN to be used by 
the onboard functions. This function is a ground function, 

- Communications Channel Power Measurement : whose aim is to estimate 
and send to the ground the total power available at output of the digitization 
chain of the communications BFN for the applied set of weighting 
operations and for the gain optimized earlier, 

- VAA processing : whose aim is to estimate the quantities J( 0 t I S( 0 t and 
J m IS u relative to the reception band B. This function is a ground function. 

For a verification by station, the method carries out the operations 
indicated by dashed lines in figure 7. From the ground and for a reception band B, it 
executes the following functions : 

- Communications Channel Power Measurement : whose aim is to estimate 
the total power available at output of the digitization chain of the 
communications BFN for the applied set of weighting operations and send 
it to the ground. This function is an onboard function parametrized from 
the ground by the function Onboard Param VAA function, 
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- VAA Gain : whose aim is the optimizing, from the results of the 
Communications Channel Power Measurement function, of the gain of the 
digitization chain of the output of the communications BFN to be used by 
the onboard functions. This function is a ground function, 

- Communications Channel Acquisition : whose aim is the acquisition and 
sending to the ground of the samples available at output of the digitization 
chain of the communications BFN for the applied set of weighting 
operations and for the previously optimized gain. This function is an 
onboard function parametrized from the ground by the Onboard Param 
VAA function, 

- VAA processing : whose aim is to estimate the quantities /„ / S u relatives to 
the stations insert. This function is a ground function. 

Onboard VAA Param Function 

Upon reception of the request for verification of the efficiency of the 
anti-jamming operation, the Onboard VAA Param function is launched. The role of 
this function is to prepare the parameters necessary for the Communications Channel 
Power Measurement functions or Communications Channel Acquisition functions. 
These parameters are : 

- the identifier of the considered coverage of the satellite, 

- the identifier of the frequency channel of the band B considered, 

- the gain of the digitization chamber of the output of the communications 
BFN to be used by the Communications Channel Power Measurement 
functions or Communications Channel Acquisition functions. Nominally, 
this gain is settled at its minimum value, 

- the function to be launched : Communications Channel Power 
Measurement functions or Communications Channel Acquisition. 

Communications Channel Power Measurement Function: 

The Communications Channel Power Measurement function has the following goals : 

- estimating the available power at output of the digitization chain of the 
communications BFN (expression (25)), 

- sending the result of the ground. 
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Communications Channel Acquisition 

The Communications Channel Acquisition function has the following aims : 

- acquiring the samples at the output of the communications BFN, 

- sending the samples to the ground. 
VAA Gain Function 

On the basis of the results of the Communications Channel Power 
Measurement function, the VAA Gain function is aimed at optimizing the gain, G x , 
of the digitization chain of the output of the communications BFN so as to exploit the 
dynamic range of encoding of the ADC to the maximum extent without saturating it. 
More specifically, this gain is computed from the results of the Communications 
Channel Power Measurement, P outpu t, of the initial gain of the reception chains, G init , 
and of the characteristics of the ADC (Gain of the ADC G adc , Maximum permissible 
power at input with the margin of 10 dB taken into account P max ). 

If the gain of the digitization chain is considered to be necessarily 
included between X and Y dB, then the function implements the following 
processing operations : 

- Computation of the power associated with the input of the ADC: P input = 

P output I Gadc , 

- Comparison of P input and P max :AP = P max - P input , 

- Computation of the gain of the digitization chain 

- If AP > 0, G x = lnf[G init + AP,Y dB] 

- If AP < 0, G x = Sup[G init + AP, X dB] 
VAA Processing Function 

For the channel verification mode, the VAA Processing Function 
implements the operations described in paragraphs IV.A to IV.F. 

For the station verification mode, the VAA Processing Function 
implements the operations described in paragraph IV. G. 



